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We investigate the dependence of the current-phase rela- 
tionship on the orientation of the order parameter for a pin- 
hole between two ^He-B reservoirs. We show that, due to 
the internal spin structure of the superfluid, the energy of 
the junction may have a relative minimum at phase difference 
equals tt at low temperatures. The dependence of the super- 
current on the direction of an applied magnetic field can be 
used to verify the present mechanism for the "7r-states" . 
PACS number: 67.40. Rp, 67.57.-z 

Superfluid '^He is a remarkable state of matter. It ex- 
hibits two superfluid phases A and B in zero magnetic 
field. In particular, in the B-phase, the magitude of the 
gap is independent of momentum direction p despite the 
fact that the pairing is triplet. It does so by having S = \ 
pairs with zero spin projections [ ^(| ti> + I iT>)] 

along direction d{p), with d related to p by a rotation: 
'^i{p) — Here R{fi,9) is the rotational ma- 

trix with rotational axis fi and angle 9. In the bulk 9 is 
determined by minimizing the dipole interaction energy 
and is given by 9^ = cos^^(— |;) « O.SStt, usually referred 
to as the Leggett angle. For any quantization axis per- 
pendicular to d, there are only | tT> and | ii> pairs. In 
the bulk the amplitudes for these pairs are equal, though 
the phase angles can be different. [Q 

Being a superfluid, one naturally expects a Josephson 
current can flow across weak-links between two reservoirs. 
This was studied experimentally flrst by Avenel and Var- 
quaux 1^ . Recently the Berkeley group have studied this 
Josephson effect again in much more detail, in particu- 
lar the current phase relationships have been mapped out 
Q . In their experiment the weak-links consist of a large 
number of small apertures (diameter « 0.1/im) made on 
a thin membrane. A particular interesting feature is the 
TT-states which occur when the temperature T is not too 
close to Tc- For an ordinary junction between two s- 
wave superconductors, the current-phase relationship is 
a slanted sine function, with the current / positive for 
phase difference x satisfying < x < tt. The junction 
energy related to the current / by / = dE/dx, is max- 
imum at X = Ti"- However, the current-phase relationships 
in H] are slanted sines only at T not too far from Tc. At 
lower temperatures, though / is positive for small x as 
usual, it turns negative at a phase difference less than tt, 
and / passes through again at % = tt with dl/dx > 0. 
Thus vr is a relative minimum instead of maximum in the 
junction free energy. Explanation of this phenomenon is 
still controversial. There is a suggestion Q that this vr 
state is not an intrinsic property of a single junction but 



rather the collective behavior of many pinholes. Earlier 
theoretical works have also predicted possible exis- 
tence of TT-states, but they rely on flnite length or the 
width of the channels. 

Later experiments reveal that there are in fact 
two possible current-phase relationships. They can be 
achieved from different cooldowns from the normal state. 
These two states are distinguishable also by the different 
magnitudes of the critical current. The authors of Ref 
1^,01 suggested that the two different states may be the 
result of two different relative orientations of fi on the 
two sides of the junction; such as parallel or anti-parallel. 

The geometry of the individual apertures in the exper- 
iment of 1^,0 approaches that of pinholes, i.e., apertures 
that have dimensions much less than the coherence length 
(w 0.1/xm). Though this criterion is not strictly obeyed 
for the experiment, pinholes are much easier to study the- 
oretically, since in this case all the self-consistent fields, 
including superfluid pairing and Fermi liquid effects, are 
the same as those near an impenetrable wall. |q| The 
current-phase relationship of a single pinhole in^He-B 
has been investigated by Kurkijarvi. He implicitly as- 
sumed that n on both sides of the junction are parallel. 
Ignoring surface depairing he found that I{x) is simply 
that of an s-wave superfluid/superconductor. We shall 
reconsider the single pinhole junction, but allowing gen- 
eral relative orientations of n'''' on the left and right of 
the junction. In particular, we shall show that one has 
a natural mechanism for the 7r-state if n' and are not 
parallel. The basic origin of the 7r-state is due to the 
internal spin structure of the order parameter. For a 
given momentum direction p, quasiparticles of different 
spin projections actually see different effective phase dif- 
ferences across the junction; thus contributing to I{x) 
with a phase shifted from each other. Provided T is not 
too close to Tc, the resultant current-phase relationship 
is anomalous, and in general can have an energy-phase 
relation that has a relative minimum at x = tt. Thus we 
have one of the very unusual situations where the spin 
structure of the order parameter affects the dynamics of 
the mass flow. 

The Berkeley experiments were performed in the 
absence of magnetic fields. In this case fi prefers to lie 
along (or opposite to) the normal near a surface thus of 
the membrane seperating the two reservoirs. However, 
under a magnetic field H its orientation can be modified. 
Ip^ We shall show that /(x) can be changed substan- 
tially by applying a magnetic field of sufficient magni- 
tude (^50G) along a general direction. Moreover, we 
predict that if one performs cool-downs from the normal 
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where Xp = X ~ fX|; ^ the area of the pinhole, Nf 
the density of states per spin at the fermi energy, Vf the 
Fermi velocity. To complete the calculation we only need 
to find xp for given n'''', with d\{p) = Rif^lji} ,6l)pii and 
similarly for I ^ r. 

In the absence of any other orientation effects such as 
magnetic field, n''*" are expected to lie along ±z. If 
n} — n"^ then obviously d} ■ (F = 1 hence Xp — ^ for 
all p. Our result for the current reduces to that of an 
s-wave superconductor ||]. For ease of later comparison 
we plot the current-phase relationship in Fig ^. This 
is the configuration with a higher critical current. Now 
consider n} = —ff. Parametrizing p by its azimuthal 
and polar angles {ap,(3p), the corresponding angles of 
d^'^ are obviously {ap ± OL.fip), One then easily gets 
d^ ■ dr ^ I - 2s\T?(3ps\T?eL = 1 - ^sin^/Jp. The resul- 
tant I{x) is as shown in Fig |^. This is the configuration 
with a lower critical current. Except for T very close to 
Tf, where I{x) is basically sinusoidal with a slight tilt, 
TT-states are evident. These current-phase relationships 
resemble closely those obtained experimentally 0| for the 
"low critical current" state. 

It is, however, known that the orientation of n can 
be affected by a magnetic field. The relevant terms in 
surface free energies are proportional to — (i • fi)"^ and 

— {HiRi^z^Y 0' "^^^ ^'"^^ term prefers n = ±i. How- 
ever, for sufficiently large magnetic field ( ^SOG) the sec- 
ond term dominates which tends to orient n in a direction 
such that the rotation Oj^) rotates ±z to H . For sim- 
plicity in the following we shall consider this case only. 
Without loss of generality we let H be in the y — z plane 
and denote its angle with the z axis hy 9h\ < Oh < 
Then the possible orientations of fi are 



state under a magnetic field in general directions, one can 
have more than two current-phase relationships possible. 
These predictions can be used to distinguish among the 
different hypotheses suggested for the 7r-state. 

We shall then consider a single pinhole between two 
reservoirs I and r with z being the direction along the in- 
terface normal. The current can be calculated along the 
same lines as in Ref. ||,^. Without loss in generality we 
take the phase of the order parameter to be and x for 
the left and right reservoirs respectively. To obtain the 
current we need to solve the Andreev |0 equation (suit- 
ably generalized to triplet pairing) or the quasiclassical 
Green's function. For simplicity we shall ignore sur- 
face depairing. Under this approximation the problem 
simplifies enormously by the observation that one only 
has S -w = ±1 pairs along any direction w perpendicular 
to d. Ti}''^ can be considered as constants in the present 
calculations since the size of the pinhole is much less than 
that of the coherence length which is in turn much less 
the bending length of the h vectors. Thus for given p and 
thus quasiparticle path through the pinhole, d is piece- 
wise constant and equals either or (F . By choosing the 
spin-quantization axis along w = d^{p) x dr{p), the gap 
matrix is finite only for the IT and ii components. Ex- 
plicitly, with the triad (m, v, w)p as the basis vectors for 
d{p), the gap matrix has the form 

A ^ f -du + idv \ . f -e-'^^ \ 
d,,+zdj = e^^" ) 

where is the azimuthal angle of d in the (u, v) plane. 

= (f>p{(f>p) for z < (>)0. The Andreev equation or the 
quasiclassical equation block-diagonalized in spin-space, 
resulting in two matrix equations only in particle-hole 
space. Each of them can be solved as in the s-wave case. 
For given p, an f quasiparticle effectively sees a phase 
TT — (pp ioi z < and ir — (p^ + x for z > 0. i.e., an 
effective phase difference of x — {4>p ~ 'l^p) ■ Similarly the 
effective phase difference for a J, quasiparticle is x + i<Pp — 
(f>p). For future convenience we shall define Xp = (f'p ^ 
(j)p. Obviously Xp corresponds to the angle between (r{p) 

and d\p), thus Xp ~ cos^^ {d^{p) ■ rf'(p)). We see that 
the contribution of the Dresent ouasiDarticle oath to the 
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With this, we can immediately see a mechanism for the 
formation of the Josephson 7r-states if xp 0, see Fig ^ 
I p^ It remains to sum over the contributions from all p. 
This can easily be done with the final result 



We shall use the letters A, B, C, D to denote the differ- 
ent orientations of fi. A and B exist only for co&Oh > ~\^ 
whereas C and D exist only for cosOh < \- At = 
the configurations A and B correspond to n = ±z respec- 
tively. A and B rotate z to H whereas C and D rotate 
—z to H. For the junction we shall denote the order 
parameter configurations on the two sides by the order 
pairs AB etc where the letters indicate n'''' respectively. 
Thus ii {) < Oh < 0.427r then the allowed configurations 
of the junction are AA, AB, BA and BB; whereas for 
0.427r < Oh < 0.527r sixteen configurations are allowed. 
The current-phase relationships of some of these config- 
urations are identical by symmetry considerations alone. 



A rotation of tt around the x axis effects the transforma- 
tions A B, C *^ D and simultaneously interchanges 
I and r. Thus e.g. Iac{x) = -Idb{-x) = -^db(x)- 
In our present approximation of no surface pair break- 
ing, I{x) depends only on d}{p) ■ dr{p). Thus we have 
^aa(x) = ^bb(x) and /ab(x) = ^ba(x) etc. We are 
thus left with 5 independent current-phase relationships 
for AA ( = BB = CC = DD) |l| AB (= BA), AC 
CA = BD = DB), AD ( = CB = BC = DA) and CD 
DC). It turns out there is also a rather non-trivial 
relation between AC and AD in that for any given p in 
AC, there exists another p' related by rotation about z 
such that d^{p) ■ dP{p) = d^{p') ■ (P{p') (Appendix A). 
Thus /ac(x) = -fAD(x)- Summarizing, for given 9h with 
Q < 9h < 0.427r there are two possible I[x)i we shall la- 
bel them AA and AB; for 0.427r < Or < O.SStt there are 
four possible /(x)'s. We denote these by AA, AB, AC 
and CD. Results for O.SStt < 6h < tt can be obtained 
from those ofO < 9h < 0.427r hy Oh ^ tt - 9h- 

As an example we show in Fig. ^ the current-phase 
relationships for these configurations at 9h — 0.457r, 
T = O.lTc. Jaa is the same as that between two s-wave 
superconductors since n'"'' are parallel. AB has energy 
minimum at x = but also a relative minimum at x ~ tt. 
AC has a rather conventional shape except for the phase 
shift by TT, thus having its energy minimum at % = vr 
rather than 0. CD has a very weak relative minimum at 

TT. 

At Oh — tt /2 the system possesses an extra symmetry: 
a rotation of tt around the z axis induces the transforma- 
tions C and B <->• D. Thus at Oh = tt/2 Jab and /cd 
merge and only three possible /(x) remain, (not shown) 

The above provides a possible test of the hypothesis 
that the tt state is the result of relative n''"" orienta- 
tions. If one performs cool down from the normal state 
in a magnetic field (of suitable orientation), in principle 
all configurations are reachable. There should be two 
possible I{x) for Q < Oh < 0.427r but at least four for 
0.427r kOh < 0.587r. [except Oh = tt/2 ] 

Next we consider the evolution of /(x) as function of 
Oh for a given configuration. We shall in particular dis- 
cuss the case where 9h is increased from 0. For AA, tt}'"^ 
remains parallel and thus /(x) is independent of 9h. The 
result for AB is as shown in Fig ||. Note as mentioned 
9h = corresponds to n antiparrallel and along ±z. As 
9h increases from initially the critical current varies in a 
non-monotonic way (Appendix B); then I{x) evolves to- 
wards the s-wave result, reaching it at 9h ~ O.SStt where 
n}'^ become parallel and both along —x. 

This provides yet another test whether the 7r-states ob- 
served in ref are due to n'"'' opposite to each other. 
Starting from the configuration where the low critical 
current state is observed, if one applies first a magnetic 
field along z of sufficient magnitude and then rotates the 
magnetic field away from the normal, provided no sud- 
den rearrangment of n'^'' takes place, the critical current 
should evolve according to Fig 0; in particular ultimately 



it should increase, the energy relative minimum at x = tt 
should become more and more shallow and eventually 
disappear near 9h = O.SStt. 

For completeness we also mention the 9h dependences 
of other configurations. /ac(x) is 9h independent under 
the present approximation (Appendix A) . /cd (x) can be 
obtained from /ab by substituting Oh ^ tt — Oh- 

Though the relative orientation between n}'"^ provides a 
natural mechanism of 7r-states, not all features observed 
in the Berkeley experiments are consistent with the 
results here. At zero magnetic field the theory here ex- 
pects TT states only for n anti-parallel (and along the nor- 
mal to the interface). We therefore must identify the 
result of Ref [|j as due to this configuration. The value 
of lo defined in the caption of Fig || corresponds to a mass 
current of ~ 2 x g / sec. Thus the critical current of 
the "low critical current state" at, e.g., T — Q.2STc is 
expected to be only around 3 x g / sec according to 
Fig^; whereas the experimental value is '--^ 7 x g / sec. 
1^ This discrepancy may be due to the finite size of the 
apertures and remains to be understood. Anyway the 
prediction of strong H dependence of /(x) here can serve 
as an important test of the hypothesis that the 7r-states 
observed are due to internal spin structure of the super- 
fluid. 

Appendix ^ - In this Appendix we discuss <f{p) ■ d^p) 
for configurations AC and AD. Obviously this dot prod- 
uct is given by PiR-ip^p^ where is the rotational ma- 
trix formed by [R{nr ,0l)]~^ R{n\OL). 7^ is thus the 
combined action of R{fi'',OL) and then the inverse of 
R{n^,OL). First we observe that since R{n\OL) rotates 
z to H whereas R{tV,Ol) rotates —z to H, TZ rotates 
z to —z. From the expressions for R{ttT'^,Ol) one can 
easily evaluate the rotational angle Q associated with Tt 
by the formula TrT^ = (1 + 2cos8). After some straight- 
forward algebra, one can obtain = tt. Thus TZ must 
correspond to a rotation of tt around an axis in the x — y 
plane. TZ for AC and AD differ only by the direction of 
this axis. Thus for any given p for AC, there exists an- 
other p' related to p by a rotation around z such that 
d^{p) ■ dP{p) = d\p') ■ dP{p') and thus their /(x) are 
identical. Also, Oh affects only the direction of the rota- 
tional axis for TZ. Thus I{x) for these configurations are 
independent of H. 

Appendix B - Here we discuss the non-monotonic 
dependence of the critical current for the configura- 
tion AB under increasing Oh. Using considerations 
along the same lines as in Appendix A, we see that 
TZ now leaves z invariant and thus TZ must correspond 
to a rotation around z itself. Q can be evaluated to 

be cos-i |i - l|. The quantity in the 

bracelets and thus is non-monotonic in Oh'- at Oh = 0, 
O = cos~^(— |); at Oh — tt/S, O has its maximum 
value of TT but then decreases upon further increase of 
Oh, reaching Q — Q &t Oh — O.SStt. Since Q is related 



to the shift of the contribution of the two different spin 
species from each other, the non-monotonic behavior of 
& results in the non-monotonic dependence of the critical 
current on 9h as shown in Fig 
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FIG. 2. Current-phase relationships for n''"" both 
along the normal and parallel to each other. The tem- 
peratures are, for decreasing critical current, T/Tc = 
0.1, 0.3, 0.5, 0.7, 0.8, 0.9. = -kANsVjI^bII 
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FIG. 3. Current-phase relationships for n'"'" both 
along the normal but opposite to each other. The tem- 
peratures are, for decreasing critical current, T/Tc = 
0.1,0.3,0.5,0.7,0.8,0.9. 
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FIG. 1. The basic mechanism of producing a 7r-state. 
Dashed and dot-dashed lines are two current-phase rela- 
tionships shifted from each other by equal and opposite 
amount horizontally. They correspond to the two terms in 
eq (|^) [shown here at T = 0]. The resultant l(x), full-line, 
is anomalous. The corresponding junction energy, being 
proportional to the integral of / over (line with dec- 
orated with symbols) has a relative minimum at x = i". 
This mechanism is operative so long as T is not too close 
to Tc, so that the individual terms in expression (^) is not 
strictly sinusoidal. 
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FIG. 4. Current-phase relationships for Qu = 0.457r. 
T = O.lTc 




A 



FIG. 5. Current-phase relationships for the AB con- 
figuration as a function of Oh/'^ given in the legend. 
T = O.lTc. 



